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Abstract 

We extend the model of ferromagnetic superexchange in dilute magnetic semiconductors to the 
ferromagnetically ordered highly insulating compounds (dilute magnetic dielectrics). The intrinsic 
ferromagnetism without free carriers is observed in oxygen-deficient films of anatase TiC>2 doped 
with transition metal impurities in cation sublattice. We suppose that ferromagnetic order arises 
due to superexchange between complexes [oxygen vacancies + magnetic impurities], which are sta- 
bilized by charge transfer from vacancies to impurities. The Hund rule controls the superexchange 
via empty vacancy related levels so that it becomes possible only for the parallel orientation of 
impurity magnetic moments. The percolation threshold for magnetic ordering is determined by the 
radius of vacancy levels, but the exchange mechanism does not require free carriers. The crucial 
role of the non-stoichiometry in formation of the ferromagnetism makes the Curie temperatures 
extremely sensitive to the methods of sample preparation. 
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I. INTRODUCTORY NOTES 



Recent experimental and theoretical investigations of dilute ferromagnetic semiconductors 
are concentrated mostly on II- VI and III-V materials doped with Mn impurities (see, e.g.,—). 
The theoretical explanation is based on the Vonsovskii-Zener model, which implies existence 
of a direct exchange interaction between the localized spins of magnetic ions (Mn) and 
itinerant spins of free carriers (holes in case of (Ga,Mn)As and related materials). 2 This 
interaction generates RKKY-type indirect exchange between Mn ions, and the latter is 
claimed to be the source of long-range FM order in dilute magnetic semiconductor (DMS) 
alloys. It was shown recently that the kinematic exchange interaction (specific version of 
superexchange between localized Mn moments via empty valence states of the host material) 
arises in p-type III-V DMS. 3 This mechanism works together with the RKKY interaction 
because of the noticeable hybridization between d-electrons of Mn ions and p-holes near the 
top of the valence band. Both mechanisms are characterized by a direct proportionality 
between the carrier concentration and Curie temperature Tq- However, it was pointed 
out recently 4 that an application of the RKKY based model to disordered DMS may be 
questionable if one properly takes oscillating character of the RKKY interaction into account. 

Meanwhile, another family of dilute ferromagnetic alloys was discovered during the recent 
five years, where the ferromagnetic order with a high Tq exists in spite of the low carrier 
concentration. The ferromagnetic order was observed even in the insulating materials with 
carriers frozen out at T — > 0. These are metal oxides doped with transition metal ions 
(see, e.g., Ref. for a review). The absence of free carriers leaves no room for a RKKY 
interaction in this case, so the question about the origin of ferromagnetism arises anew. 

The metal oxides like ZnO, SnC>2, TiC>2 are classified as wide-gap semiconductors. In 
this sense they are close to the wide-gap III-V nitrides GaN and A1N. We analyzed the 
case of n-(Ga,Mn)N in Ref. y, where the leading interaction mechanism is the Zener type 
double exchange^ vie empty states in the impurity band. Magnetic ions (MI) in all these 
compounds substitute for metallic cations. From this point of view, the difference between 
the two groups is in the charge state of substitutional impurity: the neutral state of an 
impurity should be MI 3+ in nitrides and MI 4+ in oxides. This difference implies different 
structures of spin multiplets, which is important for the magnetic properties of metallic 
alloys. However, due to a noticeable covalency of the host materials, MI should saturate the 
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broken bonds by its own valent electrons, which are donated by the 3d shells. This means 
that the violation of the electronic structure inserted by the dopants cannot be ignored 
in the studies of magnetic properties of DMS. Moreover the recent paper— emphasizes the 
role of bound exciton states in ZnO bases compounds and indicates the correlation between 
between the experimentally observed chemical trends in ferromagnetism and the electron 
binding energies to various MI. 

Among the most salient features of magnetism in dilute ferromagnetic oxides one 
should mention an extreme sensitivity of the magnetic order to the growth and anneal- 
ing conditions.-^ We believe that this is an integral feature of magnetism in these materials, 
and the ferromagnetic ordering with high Tq is mediated by intrinsic or extrinsic defects, 
which form complexes with magnetic dopants. Such point of view is supported by recent 
experimental studies of Co- dopeo- and Cr-dopecUS Ti02- It was noticed, in particular, 
that in the most perfect (Ti, Cr)02 crystals the hysteresis loop at a given temperature is 
essentially less distinct than in 'bad quality' samples.— Besides, the enormous scatter of 
effective magnetic moment per cation is observed in all ferromagnetic oxides, and the con- 
centration of magnetic dopants lies far below the percolation threshold x c associated with the 
nearest-neighbor cation coupling. 8 Basing on these facts, one concludes that non-magnetic 
defects are also involved in formation of the localized magnetic moments. They influence 
their magnitude and localization, and this influence is sensitive to the growth and annealing 
regimes. 

Recent theories of magnetism in n-type dilute magnetic oxides appeal to magnetic po- 
larons as mediators of indirect exchange between magnetic dopants.— 1 ^ In this case the 
extrinsic defects are donor impurities, which donate free electrons occupying the bottom of 
conduction band. These electrons form shallow spin-polarized polaronic states due to strong 
exchange with the localized magnetic moments of transition metal dopants. A large radius 
of these states makes the polaronic percolation threshold essentially lower than x c . 

One should note, however, that the carrier concentration in these compounds is vanish- 
ingly small, and that is why they have been proposed to be called dilute magnetic dielectrics 
(DMD). 9 It means that we need now to find an exchange mechanism, which works in the 
case of really insulating oxides. In this paper we construct a model of insulating oxide, 
which considers intrinsic non-magnetic defects as an integral part of the ordering mecha- 
nism. We concentrate mainly on the superexchange between magnetic ions mediated by 



3 



oxygen vacancies, which is believed to be the source of ferromagnetic order in Ti02 diluted 
with Co.— 

II. MODEL OF DILUTE MAGNETIC OXIDE 

As has been established three decades ago , 12 i 13 the basic microscopic Hamiltonian, which 
gives an adequate quantitative and qualitative description of electronic and magnetic prop- 
erties of semiconductors doped by transition metal ions is the Anderson Hamiltonian^ (see 
summarizing monograph o 15 ' 16 ) . Indirect RKKY 17 , superexchange and double exchange 3 
interactions between magnetic ions may be derived microscopically from the multisite gen- 
eralization of the Anderson model.— The Hamiltonian of this model reads 

H = H h + H d + H M , (1) 

where H/j describes the electronic structure of the host semiconductor, = H<# is the 
Hamiltonian of strongly interacting electrons in the (i-shells of magnetic ions located in the 
lattice sites j, and Hhd = 52,- ^hdj is the hybridization Hamiltonian describing the cova- 
lent bonding between the impurity 3c?-electrons and the host p electrons (the contribution 
of hybridization with the host s electrons in the effective exchange is generally negligibly 
small). Usually the direct dp exchange is neglected in the Anderson model, since the indirect 
exchange coupling dominates in the covalent materials . 15 > 16 

From the point of view of the Anderson model, the crucial difference between the narrow- 
gap and wide-gap semiconductors is in the position of d-levels of the impurity electrons 
relative to the top of the valence band. Due to a strong electron-electron interaction in the 
3d shell (Coulomb blockade), the impurity energy level Ed should be defined as an "addition 
energy" for the reaction d n ~ x + e& — > d n , i.e. the energy cost of adding one more electron, 
e&, from the band continuum to the impurity 3d-shell, 

e d = E(d n ) - E^- 1 ) (2) 

where E(d p ) is the total energy of a "pseudo-ion" consisting of an impurity in a configuration 
3d p and distorted electron distribution in the host semiconductor 15 . 

In this definition the energy is counted from the top of the valence band, which is the 
source of additional electrons for the impurity. It is known, e.g. that Mn, the most popular 
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magnetic impurity, has especially deep d level corresponding to the half-filled d shell in 
the configuration 3d 5 . This level is always occupied in such narrow-band semiconductors as 
(Ga,Mn)As, (Ga,Mn)P, (In,Mn)As. This means that the addition energy e d = E(d 5 )-E(d' lt ) 
is well below the energy, e p , of the top of the valence band. On the contrary, > e p in the 
wide gap semiconductors like GaN,— and the occupation of this state depends on the type 
and concentration of other defects. 

One should compare this picture with the energy spectra of magnetically doped metal 
oxides. Situation with the transition metal impurities in ZnO is apparently close to that 
in GaN. Both these materials are wide band gap semiconductors, and the only difference 
is that the neutral state of magnetic impurity in ZnO is MI 2+ . The corresponding energy 
levels rule, within the band gap, and the states with other oxidation numbers 

may appear only due to an additional n - or p - doping. 

The calculation^ of addition energies for MI n+ states in rutile Ti02 has shown that the 
"energy levels" Ed for neutral MI 4+ states are usually deep in the valence band for heavy 
transition metal ions (starting with Mn), whereas the condition Ed > s p is fulfilled for 
the addition energies of charged MI 3+ and MI 2+ ions. Recent LSDA+U calculations^ for 
anatase Ti02 doped with Mn, Fe and Ni are in a general agreement with those earlier cluster 
calculations. 

This means that the ground state of this system should be Ti02:MI 4+ . In this situation 
the Anderson Hamiltonian may be immediately transformed into the sci-exchange model by 
means of the Schrieffer- Wolff transformation. In the absence of free carriers the short-range 
effective inter-impurity exchange interaction is exponentially weak and there is no chance for 
magnetic ordering with reasonably high Tc- However, the situation in real materials is quite 
different. Leaving aside early results for the samples with secondary phase inclusions ,— we 
concentrate on the recent data for the materials, which are believed to be genuine ferromag- 
netic DMD free from magnetic precipitates. As was mentioned above, the room temperature 
ferromagnetic order was detected in Cr and Co doped anatase Ti02, and in both cases na- 
tive defects are involved in formation of the ferromagnetic order.^i In the latter case these 
defects are apparently oxygen vacancies. Below, we adapt our theory^ for vacancy mediated 
indirect exchange. The case of (Ti,Cr)02 will be discussed in the concluding section. 
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III. VACANCY-RELATED SUPEREXCHANGE MECHANISM 



In the highly insulating ferromagnetic (Ti,Co)02, magnetic ions are in Co 2+ state.— A 
seeming contradiction with the requirement of electrical neutrality of the substitutional 
impurities 2 ^ may be resolved, if one takes into account the intrinsic non-stoichiometry of the 
samples: oxygen vacancies Vo, are created in the process of sample preparation. It follows 
from the general neutrality consideration that an oxygen vacancy binds two electrons on a 
discrete level in the band gap to saturate the dangling bonds. In the process of annealing a 
noticeable amount of oxygen vacancies are captured by Co substitution impurities in cation 
sites. In the case, when Vo arises in the octahedron containing a Co ion, the double defect 
Co — Vo is formed, and its stability (or metastability with respect to formation of Ti — Vo 
complexes) is determined by the charge transfer of these two electrons in accordance with 
the reaction 

V (p 2 ) + Co 4+ (d 5 ) - [V + Co 2 +(d 7 )]. (3) 

A presumable scheme of the energy levels, which substantiates such charge transfer is shown 
in Fig. 1. The configurations d n of 3d shells corresponding to different charge states of Co 
and Cr impurities are pointed out. The dashed band shows the position of vacancy- related 
levels below tho bottom of conduction band e c . 

As a result of charge transfer reaction (JHJ), the vacancy related electrons are captured by 
[V 2+ Co 2+ ] compl exes, the system remains insulating and the magnetic interaction between 
these complexes predetermines the magnetic properties of the system. 

In this section we construct a microscopic model of impurity-vacancy complexes basing 
on our previous studies^ of double ferromagnetic exchange in (Ga,Mn)As and the electronic 
structure of complexes Cr-Vo in (Ga,Cr)As,— as well as on the numerical calculations of the 
electronic structure of Vo defects in TiOjrS and doped non-stoichiometric (oxygen deficient) 
(Ti, Co)0 2 ^. 

The structure of Hamiltonian, which describes the oxygen deficient (Ti,Co)02 is the 
same as in Eq. (fl|). but the oxygen vacancies should be taken into account in the host 
Hamiltonian Hh- We start with a quantum-chemical picture of an isolated defect cell, which 
contains both Co impurity substituting for a Ti cation and the vacancy in one of the apices 
of oxygen octahedra. We mark the cell by the index j and consider the Hamiltonian 

Hvj = W h j + + Hhdj- (4) 
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FIG. 1: Impurity-related and vacancy-related energy levels in cation substituted non-stoichiometric 
Ti0 2 . 

Now the first term )M h - describes the host crystal with the vacancy, Vo, in the cell j. This 
Hamiltonian may be diagonalized, and the resulting spectrum contains the band continuum 
distorted by the vacancy potential and the discrete levels e Vfi for electrons bound on the 
dangling bonds in the defect cell. 

The Hamiltonian rewritten in this basis contains terms describing hybridization with 
both continuum and discrete states of H^ -. The general form of H hd j is 

H hdj = ^^(Vju^T d ,T )d] +Td41(7 c j+To , UCT + h.c^ 

\wa T d T 

+ Yl Yl (y^ Trf ' K H+T d ,»* c ™° + h - c ) ■ ( 5 ) 

k/^<x a=v,c 

Here the operators dj +Td ^ a and Cj +TotUa stand for the localized d\x - and pu - orbitals on 
the impurity site Rj = Rj + and the vacancy cite R D = Hj + r Q , respectively, and 
the operators c Kaa describe the continuous states in the host crystal valence (a = v) and 
conduction (a = c) bands. 

Information about the vacancy-related states in anatase TiC>2 is rather scanty, 28 so we 



refer to the data available for oxygen-deficient rutile Ti02- In accordance with numerical 
calculation o 29 i 30 , which correlate with the experimental observations, Vo creates donor levels 
e a few tenth of an eV below the conduction band. The donor electrons saturate the dangling 
bonds with the neighboring atoms, so that the defect wave function has the largest amplitude 
at the next Ti neighbor of the vacancy and extends to several coordination spheres. The 
extended defect states form a band of donor levels already for 1% concentration of oxygen 
vacancies, i.e. the wave functions of the vacancy states effectively overlap when the vacancies 
are separated by a distance of about 5 interatomic spacings. This observation will be used 
below when estimating the Curie temperature. The calculated electronic structure for Co- 
doped oxygen deficient anatase Ti0 2 gives similar picture for the charge distribution and 
the density of states.— 

This observation verifies the model Hamiltonian for the Co substitutional impurities. 
Its form implies a strong nearest-neighbor hybridization between the vacancy and impurity 
related states. Since the conduction band of Ti02 is formed mainly by the d-states of Ti 
sublattice, the matrix elements V CVL in the second term in this Hamiltonian describe the 
hybridization of the Co impurity rf-states with those in the cationic sublattice. Apparently 
one may neglect the hybridization V Cfl in comparison with Vj ulM , and retain in this term only 
the hybridization Vj Vlx with the oxygen-related p-states in the valence band. 

We return now to are the energy level scheme Fig. ^substantiating the charge transfer 
reaction (J3J). It is emphasized that we compare the one-electron energy levels of the oxygen 
vacancy with the corresponding levels of Co impurities in different charge states. These levels 
are defined as addition energies (J2J) corresponding to recharging processes Co m+ /Co <m ~ 1 ^ + . 
These energies may be taken, i.g. from Ref. |20|. Hybridization Vj Vfl allows for an electron 
transfer of the two electrons from the vacancy level first to the level Sd{3 + /4+) and then to 
the level £^(2 + /3+) and formation of the complex defect [V G + Co 2+ ]j with seven electrons 
occupying the 3d-shell of the ion Co 2+ and empty extended Vo-related states. 

Indirect exchange interaction arises when the wave functions of two complexes centered 
in the cells j and i overlap. It is clear that the overlap is controlled by the radius of 
extended Vo-related state. An obvious interaction mechanism is the superexchange, where 
two magnetic impurities exchange electrons via empty vacancy levels. It will be shown below 
that this interaction favors ferromagnetic ordering in the system. 

The simplest way to construct an effective superexchange interaction operator is to 
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project the two-impurity Hamiltonian Hjj = + Hdi + H' hji + Hhdj + ^hdi onto the subspace 
. . . with fixed Co 2+ configurations of the two magnetic ions and also project out 

the charged impurity states Co 3+ , which arise due to the hybridization H hd j + H^. Here 
H'^j is the host Hamiltonian with two oxygen vacancies in the cells Then the fourth 
order term in the impurity-vacancy hybridization 

(ji I ^hdj p no - Hhdi p rw^hdi p no" Hhdj \]i)2+ (6) 

Eiq — r\j { Ujq — hiQ — r\ j{ 

gives the effective exchange vertex, represented graphically in Fig. |21 The resolvents, which 
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FIG. 2: Diagrammatic representation of the superexchange interaction between two magnetic 
impurities in charge states df and dj via the intermediate states d™ -1 and d™" 1 and empty levels 
pa (see text for further explanation). 



appear in (JBJ) correspond to the ground state Eq of the non-hybridized Hamiltonian = 
H<# + H di + H^j. The solid lines in this diagram stand for the definite charge and spin 
state of the impurity in the site i or j (Co 2+ (d 7 ) and (Co 3+ (<i 6 ) in a given spin and orbital 
configuration). The vertices correspond to various terms in the hybridization Hamiltonian 
(0), which changes the impurity state from dj to d® and v. v. In the vertex V v ^ an electron 
appears on the empty vacancy v - orbital with the spin a. In the vertex V afJi a hole (electron) 
appears in the valence (conduction) band. Yet, a similar process dj — > d® + v'o occurs 
around the counterpart site j. Then the two p-electrons change their 'hosts' (two dashed 
lines in Fig. EJ), and this is the end of the superexchange process. This scheme allows one 
to roughly estimate the effective exchange coupling as 

J Z ~ V*/AlD , (7) 

where A^ is the energy of the transition (J3j), and D a is the characteristic energy scale of 
the vacancy related band under the bottom of conduction band. 
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A similar estimate may be found, e.g. in Refs. j23j, where the superexchange between 
singly occupied states in the 4 /-shell via empty states in conduction band was studied (D 
in that case is the conduction band width). In the case of a singly occupied shell f 1 or 
d 1 there is no Hund interaction in the unfilled shell. As a result antiferromagnetic part of 
this interaction prevails (ferromagnetic channels are blocked by the Pauli principle).— In 
our case the superexchange mechanism favors the ferromagnetic alignment of spins since the 
intrashell Hund interaction suppresses the exchange by the electrons with antiparallel spins 
(see next section). 

There seems to be much in common between the vacancy-related mechanism represented 
by Eqs. ((7j) and the spin-polaron mechanism offered in Ref. 0. In both cases the effective 
radius of magnetic interaction is controlled by the extended donor states. It should be 
emphasized, however, that unlike the indirect exchange via the states in the donor impurity 
band&, in our case the band of vacancy states may be empty (provided all the electrons 
from the vacancy levels are captured by acceptor-like magnetic impurities), so that this 
mechanism may be realized in highly insulating materials. On the other hand, the empty 
donor band of extended vacancy related states plays, in oxygen deficient (Ti)02, the role 
similar to that played by the empty states in the hole pockets of p-type (Ga,Mn)As in 
superexchange between magnetic ions (see Ref. 0). We will return to the comparison of 
superexchange mechanism in various dilute magnetic compounds in the concluding section. 



IV. MICROSCOPIC THEORY OF SUPEREXCHANGE IN DMD 

The simple estimate (J7J) by means of the fourth order perturbation equation (JHJ) may 
be essentially improved. Basing on the theory of isolated transition metal impurities in 
semiconductors , 15 ^ 6 one may consider the effects of single-site hybridization between the 
o?-orbitals of impurity and oxygen-related states in the valence band and vacancy levels and 
then use the perturbation theory approach for the intersite hybridization effects responsible 
for the superexchange. 

It is known that the main result of this hybridization is 'swelling' of the wave functions 
of the impurity ci-electrons . 15116 These functions acquire 'tails' formed by a superposition of 
Bloch waves from valence and conduction bands. A localized vacancy state is also a super- 
position of the same Bloch waves. All these covalent effects modify the simple estimate 
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which necessitate a more refined analysis within the microscopic Anderson-like Hamiltonian. 
In this section we will derive equations for the superexchange interaction in the system 



For the 



of defects [Coxi, Vo] by means of the Green function method developed in Ref. 
purpose of practical calculations within this method, it is more convenient to write the term 
H' h j in the model Hamiltonian (jlj) in the form 

^hj = ^ ] £ nC Kaa C Kaa + ^ ] £ w C -j+T v ,ucT C j+Tv,V(T- (8) 
nacr va 

It is assumed in Eq. (jHJ) that the host Hamiltonian with vacancies has been diagonalized, 
and the band states are classified by the index k. Each vacancy in the site j + t v creates a 
localized level in the gap. It was found in the previous numerical calculation o 20 ' 29 that both 
the oxygen-related p-states in the valence band and the titanium-related d-states from the 
valence band contribute to formation of localized vacancy states, and the vacancy-related 
level arises slightly below the bottom of conduction band. 

The term describing the subsystem of magnetic impurities in the Hamiltonian has 
the form 

h, = E H <* = E E £ m4 a* + H r r - ( 9 ) 

3 3 A" 7 

Here the site indices are introduced in accordance with the definition given in (|5|). Hf orr 
includes all the Coulomb and exchange interactions responsible for the Hund rule within the 
3<i-shells of Co ions. In accordance with this rule the Co d-shell contains seven electrons, 
five of which form the closed 'inert' d^(t\e 2 ) subshell and the remaining two electrons form 
the open 'active' d^At^e^) subshell. Only the electrons in the open subshell are involved in 
the hybridization induced superexchange. The reason for this discrimination is the above 
mentioned swelling of the impurity electron wave functions due to the hybridization with 
the host electrons, which is stronger for the open subshell than for the closed one . 15 ' 16 
It is easily seen from © that the fourth order transition 

for ferromagnetically aligned spins of two impurities is possible, whereas the similar process 

is suppressed due to the Hund rule for the occupation of subshells. Here Rji = (E — H^)" 1 
is the zeroth order resolvent. 
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Our task is to calculate the energy gain due to the electron exchange between the impu- 
rities with parallel spin alignment. Since no spin-flip occurs in the course of the exchange, 
the spin index may be omitted in the Hamiltonian (JHJ), (jHJ), ©• This energy gain is given 
by the following equation*^ 

1 f°° 

-Im / £TrAG[e — i8 sign(e — rj)]de, (10) 

T J-oc 
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where AG is the part of the single electron Green function related to the spin-dependent 
interaction of magnetic impurities, 

AG(z) = (z - H)- 1 - (z - H' h] y\ 

r) is the chemical potential. The integral in (|10j) contains contributions from both the valence 
band continuum and occupied discrete states in the energy gap. 

When calculating AG(z), we exploited the following features of our problem. 

(i) Since we are interested in the energy gain due to the parallel spin ordering, the spin-flip 
processes do not contribute to the relevant part of the system of Dyson equations for the 
Green function and these equations allow for an exact solution.— 

(ii) The calculation procedure may be radically simplified by using the analytical properties 
of TrAG(^), which allow one to calculate the energy by summing over the empty states 
instead of integrating over the whole occupied part of the energy spectrum.— 

(iii) We are interested only in the trace of the Green function defined in the subspace 

(441 •••144) ( seeE q- ©)• 

As was discussed above, the leading contribution to the formation of complex defect is 

due to the hybridization between the d-shells of magnetic impurities and the empty states 

in the band of vacancy related states below the bottom of conduction band. Thus, our next 

task is to calculate the upward shift of the Vo-related levels. For this sake one has to find 

the poles of the single electron Green function G. This function has the following block 

structure 



Gfj, v ' » 
G VAt ' G vv ' 

y Gy^ Gyyl j 



(11) 



where the diagonal blocks are the two impurity (G w /), band (G w ') and single- vacancy {G w >) 
Green functions, respectively, and the off-diagonal blocks stand for two types of hybridiza- 
tion. Only valence band states are taken into account in the band component. 
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As shown above the exchange energy E magn is given by the shifts of the Vo levels which 
are determined by the poles of G. It is sufficient to study the block G uu only since it has the 
same poles. Moreover, as follows from Eq. ()1U|). we need only its diagonal part G vv in order 
to calculate its trace. It reads 

G w = g v + gl V: [i (R i )G ifliifl V u>ft (R i )+ 

i=l,2; fx 

9l{e) E K;(Ri)G^( £ )K<^ (12) 

jj'=l,2; fifi' , v' 

We have kept the leading terms in the hybridization of the vacancy and impurity states 
contributing to the intersite exchange. Hybridization of the empty conduction band states 
and the impurity states is neglected in calculation of this exchange (see discussion after Eq. 
(jSJ) in Section III). Since the spin-flip processes are not involved in the calculation of E ma9n , 
the equations may be solved for each spin projection separately, and we omit below the spin 
index for the sake of brevity. 

We represent the impurity Green function in terms of the irreducible representation of 
the symmetry group of axial defect [Coti,Vo], so that the impurity Green functions are 
diagonal in this representation 

Gd-i^ — [s — Sdyi/j, — 1 (13) 

with the mass operator 

^d;in(^) = E K*M(^-i)^( £: )K At (H-i)G'd;t,^- 

which results from the exact solution of corresponding single impurity problem for a given 
spin . 13 i 26 Here the index k = {v, z/} unites both band and vacancy indices. The bare propa- 
gators are defined as 

1 1 1 

9v j 9v j gd;ia 

£ S v E E v E E{ a 

Using Eqs. ()12|) and ()10j) we can find the change of the total energy due to the interaction 
of the vacancy and impurity states. The second order correction in the vacancy - impurity 
hybridization V^^(Rj) in Eq. (fT2|) is a single impurity effect and has nothing to do with the 
exchange interaction. The latter appears only in the last term of Eq. (|12jl . which is of the 
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fourth order and describes a shift of the vacancy levels due to the electron exchange between 
two impurities, which is possible only provided the spins of two impurities are parallel. 

The upward shift of the empty vacancy levels corresponds to lowering of the total energy 
of the system (see discussion after Eq. ()1U|)). Hence, we come to the equation for the energy 
gain due to an indirect exchange (superexchange) between the two magnetic impurities, 

E ma gn = 2 y VuWVM^WA^)Vv,A^) (14) 

ij v £/ m , [e« - - Pdyi^u)] {su - e v t)[e' v - e d .^> - P d , m , {e' v )\ 

where Pd,ia( £ ) — ReS d j a (£:). Eq. (fTlj) has a structure similar to that of eq. (10) of our 
paper 3 . The role of hole pockets is now played by the empty vacancy states, and the fact that 
there is no free carriers and, hence no Fermi surface, so important for the RKKY model, play 
now no role whatsoever. Besides, contrary to the polaronic model^, the ferromagnetism with 
high enough Tq arises in our model without additional enhancement due to spin polarization 
of the vacancy related band (see more detailed discussion below). 

The fourth order exchange energy can be interpreted by means of the diagram presented 
in Fig. El Now the dashed lines in this diagram correspond to the two vacancy states 
virtually occupied in the superexchange act. Two multipliers in the denominator in (j!4)l are 
the energies A d0t j required for the two "reactions" Co(<i 7 ) — > Co(<i 6 ) + Vo(p) in the sites 
i,j. These energies, include also the ligand field shifts Pd-i^i of the impurity d levels . 12 i 13 i 25 
The vacancy band is empty so that the summation in Eq. ()12J) is carried out over all the 
states of the band. As a result the third multiplier, (e v — e v >), in the denominator is of the 
order of the vacancy band width D Q . Thus the calculation of exchange energy by the Green 
function method by means of Eqs. (|l()jl . (fT4*|) confirms the qualitative estimate (jJJ). 

Now we are in a position to derive the effective spin Hamiltonian for interacting magnetic 
defects. Since we integrated out all the charge degrees of freedom when calculating the energy 
j^magn ^ ^ g jj airiiH.oiii an has the simple form 

(ij) 

where Sj is the spin of the ion Co(d 7 ), and the exchange constant Juj\ is determined as the 
energy difference per bond between the parallel and antiparallel orientations of spins in the 
pair (ij). One may identify the coupling constants as 2J^j~ ) S(S + 1) = E™ a9n . Since these 
coupling constants are negative, they lead to a ferromagnetic ordering at temperatures below 
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the Curie temperature Tc- Well beyond the percolation threshold (which may be identified 
with the minimal vacancy concentration sufficient for formation of the vacancy band) Tc 
may be estimated as 

Here \E magn \, and respectively | J\, is the typical value of the superexchange interaction in the 
pair of impurities separated by a distance not exceeding the double localization radius Ro of 
the vacancy state, z is the number of Co impurities within the corresponding volume. As was 
mentioned above the fact that 1% of vacancies is sufficient for formation of a band, meaning 
that the typical distance between them is smaller than 2Ro- Since nearly each Co atom forms 
a complex with an oxygen vacancy, a 2% concentration of Co atoms substitutions means 
a 1% vacancy concentration when the vacancy states form a band. This fact guaranties 
that each Co atom interacts not only with its own vacancy but also with other neighboring 
vacancies, which form pairs with other Co atoms. 

Magnetic ordering with a high Tc in (Ti, Co)02 is observed at higher Co concentrations 
up to 6% or 3% of the concentration of vacancies captured by the Co impurities. A 1% 
concentration of vacancies creates a band in the forbidden energy gap so that we may expect 
that the coordination number for such a concentration may be at least 4. At a three times 
higher concentration the coordination number will certainly become higher and a rough 
estimate z « 10 seems to be quite reasonable. We do not dispose precise experimental 
data concerning positions of the Vo level in anatase Ti02, so we refer the data available 
for the rutile modification where the donor levels related to an isolated Vo are found at a 
few tenth of 1 eV below the bottom of conduction band.— Using for the forbidden band 
width the experimental value 3.2 eV and basing on the calculations of Co-related d-level 
positions in the forbidden energy gap presented in Ref. |2l|, one may estimate the pd charge 
transfer gap A p d ~ leV. The width of the vacancy band D as well as the magnitude 
of hybridization parameter strongly depend on the specific characteristics of the sample 
(vacancy concentration, annealing regime etc). Besides, the vacancy band may or may be 
not partially filled due to uncontrollable donor impurities. Taking for an estimate the average 
value of V Vfi ~ 0.1 eV, assuming the same estimate for D Q , and substituting the value of 
S = 3/2 into Eq. (jl(i)l one gets Tc ~ jz ■ 10 _4 eV ~ 7IOOK for z = 10. The coefficient 7 
includes all the uncertainties (degeneracy factor of vacancy and impurity levels, difference 
between the intracell and intercell hybridization parameters, uncertainties in the effective 
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coordination number z etc. We expect the factor 7 to be ideally close to one. However, 
it may deviate, and even strongly, from this value depending on the preparation history of 
each particular sample. Unfortunately, it is rather difficult to find optimal conditions for 
high Tq in such a multifactor situation. 



V. CONCLUDING REMARKS 



Among many dilute magnetic dielectrics only one example has been chosen, namely 
oxygen deficient (Ti, Co)C>2, and it has been shown that the imperfection of this crystal is 
crucially important for the formation of a long-range magnetic order. This idea is supported 
by the recent experimental observation: ferromagnetism is suppressed in Co doped Ti02 
film with high structural quality— Similar correlations between the quality of films and the 
magnetism were discovered in Ref. [l0| for (Ti, Cr)02- It follows directly from our theory that 
without additional defects creating shallow levels under the bottom of conduction band, the 
radius of magnetic correlation is too short to overcome the percolation threshold, so that 
improvement of the film quality is detrimental for magnetic ordering in DMD. 



It is also noticed in Ref. 



lOl that the charge state of chromium impurities in Ti02 is Cr H 



This fact is easily understood within the framework of our hypothesis about the origin of 
magnetic order in DMD. Indeed, the difference between Co and Cr from the point of view 
of their electronic levels in host Ti02 [in fact, the addition energies defined in Eq. (J2J) for 
the two impurities] is the position of these levels relative to the shallow defect levels. In 
correlation with the experimental data on ESR signal and theoretical calculations of the 
addition energies,— the level of Cr 2+ state in rutile Ti02 is higher than the levels of shallow 
impurities and apparently higher than the bottom of conduction band. We assume that the 
situation is similar in anatase Ti02, so that the same charge transfer mechanism, which binds 
the ions Co 2+ with oxygen vacancies, should stabilize chromium impurities in a state Cr 3+ 
(see Fig. 1). It is worth also mentioning that generally we may expect that the chemical 
trends in the Curie temperature for various transition metal impurities should correlate with 
the deep level energies (e.g.,—). 

Another experimental fact, which demands theoretical explanation is a strong sample 
dependent scatter of saturation magnetic moment M s , which also depends on the fabrication 
conditions4*i£ This fact may be easily explained if one recognizes that only part of oxygen 
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vacancies is bound to magnetic impurities. "Free" vacancies, Vo, retain their electrons and 
eventually donate them to the vacancy band. The partially occupied band is spin polarized, 
and this polarization may enhance or partially compensate the magnetization of transition 
metal ions depending on the net sign of effective exchange in this band . 8 ' 28 ' 30 

To conclude, we have proposed in this paper a mechanism of ferromagnetic exchange in 
dielectric non-stoichiometric TiC>2 doped with transition metal impurities, which involves 
the interaction between magnetic ions and oxygen vacancies and proposed a qualitative ex- 
planation of the basic properties of these materials. A quantitative calculations of electronic 
structure of magnetically doped TiC>2 based on the model described in this paper will be 
presented in forthcoming publications. 
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